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1 Introduction 

Recent high-precision measurements of forbidden transitions in highly charged argon ions [1] 
provide very good possibilities for tests of relativistic and quantum electrodynamic (QED) 
effects in middle-Z few-electron systems. In particular, in the case of Ar^^"*" the QED contri- 
bution to the 2s^ 2p — ^ -P3/2 transition energy is by four orders of magnitude larger than 
the corresponding experimental error. From the theoretical side, however, large efforts must 
be undertaken to achieve this accuracy (for the current theoretical results, see [1]). The main 
goal of this paper is to calculate the isotope shifts of forbiden transitions in Be-like and B-like 
argon and to examine a possibility for their experimental determination. 

The isotope shifts in the atomic transition energies arise from the finite nuclear mass (recoil 
effect) and from the finite nuclear size (field or volume shift). Generally, the field shift, caused 
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by the penetration of the electron wave functions into the nuclear region, dominates for heavy 
atoms whereas the mass shift dominates for light atoms. 

The relativistic units {h — c— 1) are used in the paper. 



2 Mass shift 

The nonrelativistic (NR) theory of the mass shift for many-electron atoms was first formulated 
by Hughes and Eckart [2] . This shift can be represented as the sum of two parts, the normal 
mass shift (NMS) and the specific mass shift (SMS). The normal mass shift operator is a 
one-particle operator. It is given by 

i 

where Pj is the momentum operator and M is the nuclear mass. The specific mass shift operator, 
which is a two-particle operator, can be written in the form 

^S^^EPrP.- (2) 

The sum of expressions (1) and (2) defines the total recoil operator in the nonrelativistic theory. 

The full relativistic theory of the nuclear recoil effect can be formulated only within quantum 
electrodynamics. Such a theory was first formulated in [3], where the complete aZ - dependence 
formulas for the recoil corrections to the atomic energy levels to first order mm/M were derived. 
Later, these formulas were rederived in simpler ways for the case of a hydrogenlike atom [4, 5] as 
well as for a general case of a many-electron atom [6] (see also [7]). Full relativistic calculations 
of the recoil effect, based on these formulas, were performed for H- and Li-like ions [8, 9, 10]. In 
[5], they were employed to calculate the recoil corrections in H-like atoms to order {a.Z)^w? / M , 
where a is the fine structure constant. As follows from these formulas, within the lowest-order 
relativistic approximation (~ v^jc?) and to first order in m/M, the recoil corrections can be 
derived by using the following recoil Hamiltonian 
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(3) 



where o: is a vector incorporating the Dirac matrices. The expectation value of Hm on the Dirac 
wave function yields the recoil correction to the atomic energy level to first order inm/M (here 
and in what follows, the Dirac wave functions are the eigenvectors of the Dirac-Coulomb-Breit 
Hamiltonian). An independent derivation of Hamiltonian (3) was presented in [11]. In [12], 
this Hamiltonian was employed to calculate the {aZYm^ / M corrections to the energy levels of 
He- and Li-likc middle-Z ions. 

According to expression (3), the lowest-order relativistic correction to the one-electron recoil 
operator is given by 

Tj(r) I ^aZ ( ioLi- rj)rA 
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The corresponding two-electron correction is 

' Q,. + '^'-^ ] . p (5) 




To the lowest order in m/M, the mass isotope shift is determined as the difference of the 
expectation values of the recoil Hamiltonian for two different isotopes: 

5E^s = mHu, - Hum = -^^K, (6) 

where 5M — Mi — M2 is the nuclear mass difference, 1-0) is the eigenvector of the Dirac- 
Coulomb-Breit Hamiltonian, and the constant K is defined by 

K/M = {il:\HM\^) . (7) 

It should be stressed here that our approach to the relativistic recoil effect differs from those 
in [13, 14, 15, 16], where the nonrelativistic recoil Hamiltonian was evaluated with the Dirac 
wave functions. Averaging the nonrelativistic recoil operator with the Dirac wave functions 
strongly overestimates the relativistic correction to the recoil effect. It is caused by the fact (see 
[3, 8]) that the relativistic contribution that results from averaging the nonrelativistic nuclear 
recoil operator with the Dirac wave functions is considerably cancelled by the contribution of 
the relativistic correction operator defined by equations (4), (5). For this reason, the relativistic 
wave functions must be employed only in calculations with the relativistic recoil Hamiltonian 
(3). 



3 Field shift 

The field isotope shift is caused by different nuclear charge distributions for different isotopes. 
For the nuclear charge distribution we used the Fermi model: 

N 

p{r,R) = — r, ry^, (8) 

1 + exp[(r — c)/a\ 

where the parameter a was chosen to be a = 2.3/4 In 3 and parameters N and c can be obtained 
from the known value of the root-mean-square (rms) nuclear charge radius it! = {r'^Y^'^ and from 
the normalization condition for p(r, R): 

J drp{r,R) = 1. (9) 

With a high accuracy, N and c can be determined by the following analytical formulas (see, 
e.g., [17]) 

3 ttVn-i 



|(^2)_7^2^2. (10) 
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The potential of an extended nucleus is given by 

r> = max(r, r'). (11) 

This potential was used in the Dirac-Coulomb-Breit Hamiltonian to obtain the relativistic wave 
functions. The field isotope shift can be determined by the formula 

SEps - {i^\J2^Vr,{u,Rm, (12) 

i 

where 

SV^ir) = VMir,R + SR) - \^7v(r,i?), (13) 

6R is the difference of the rms nuclear charge radii for two isotopes, and d{r^) is the related 
difference of the mean-square nuclear radii. Within the precision required, the direct evaluation 
of the field isotope shift by solving the Dirac-Coulomb-Breit equation for two different isotopes 
yields the same results. 

Concluding this section, wc note that in the case of a hydrogcnlikc atom the isotope shift 
can easily be calculated according to analytical formulas presented in [17]. 
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VN{r,R) = -inaZ [ dr' p{r' , R) — , 



4 Method of calculation 

In this work, the large-scale configuration-interaction (CI) Dirac-Fock (DF) method was used to 
solve the Dirac-Coulomb-Breit equation and to calculate the energies and the isotope shifts of 
the forbidden transitions in Ar^^+, Ar^^"*", and Kr^^+. The many-electron wave function ^(7 J) 
with quantum numbers 7 and J was expanded in terms of a large number of the configuration 
state functions (CFSs) with the same J 

*(7J) = J2^-^-{J). (14) 

a 

For every relativistic atomic configuration the CFSs $a(J) are eigenfunctions of the square of 
total angular momentum and they can be obtained as the linear combinations of the Slater 
determinants corresponding to this configuration. The set of the CFSs in expansion (14) was 
generated including all single and double excitations and some part of triple excitations. 

The Slater determinants are constructed from one-electron four-component Dirac spinors 
(orbitals). For the occupied shells these orbitals {ipj) were obtained by the multiconfiguration 
Dirac-Fock (DF) method. The other vacant orbitals {(pj) were obtained by solving the Dirac- 
Fock-Sturm equations 

<fj = ^jW{r)^j, (15) 

where h^^ is the Dirac-Fock operator, Sj^ is the one-electron energy of the occupied DF orbital 
ifjQ, and W{r) is a constant sign weight function. The parameter Xj in equation (15) can be 
considered as an eigenvalue of the Sturmian operator. If W{r) — > at r — > 00, all Sturmian 
functions cpj have the same asymptotics at r — > 00. It is clear that for Xj = the Sturmian 
function coincides with the reference DF orbital ipj^. The widely known choice of the weight 
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function is W{r) = 1/r, which leads to the well known 'charge quantization'. In the relativistic 
case this choice of the weight function is not very successful. In our calculations we used the 
following weight function 

l-exp(-(arff 

[ary 

In contrast to 1/r, this weight function is regular at origin. It is well-known that the Sturmian 
operator is Hermitian and it does not contain continuum spectra in contrast to the Fock oper- 
ator. Therefore, the set of the Sturmian eigenfunctions forms the discrete and complete basis 
set of one-electron wave functions. 



5 Results and discussion 

Before discussing our results for the isotope shifts in argon ions, we examine our calculations 
of the specific mass shift for the 3pi/2 — 3s transition in Na-like ions and compare them with 
the related results from [15]. In Tabic 1 we present our results for the SMS contribution to the 
isotope-shift constant K, defined by equation (7). The second column shows the nonrelativistic 
values of this contribution. These values were obtained within the same CIDF method and 
the same computer code, as described above, by the 1000-times increase of the velocity of hght 
(in atomic units) . We verified this nonrelativistic limit by comparing the total energies for the 
Dirac-Coulomb-Breit Hamiltonian with the total energies obtained by the fully nonrelativistic 
method based on the Schrodinger Hamiltonian and on the same calculation scheme. The total 
energies exactly coincide with each other for all the ions. In the third column we present 
the results obtained by averaging the nonrelativistic recoil operator (2) on the Dirac wave 
functions. The fourth column indicates the relativistic correction that results from using the 
Dirac wave functions in the calcTilation of the nonrclativistic-recoil-operator contribution. In the 
fifth column we give the relativistic correction obtained by averaging the relativistic correction 
operator (5). In the sixth column the total relativistic data are listed. The last column indicates 
the results of [15], based on averaging the nonrelativistic recoil operator with the Dirac wave 
functions. Comparing the relativistic corrections, presented in the fourth and fifth columns, we 
observe that the absolute value of the second relativistic correction is larger than the absolute 
value of the first one and, as a result, the total relativistic correction and the correction that 
accounts only for the relativistic effects on the wave functions are of opposite sign. This means 
that using only the nonrelativistic part of the recoil operator in calculations based on employing 
the Dirac wave functions gives strongly incorrect results for the relativistic recoil effect. 

For the 3pi/2 — 3s transition in neutral Na, the experimental value of the isotope shift 
^^Na/^^Na amounts to -758.5(0.7) MHz [18]. To derive the corresponding SMS we have to sub- 
tract the NMS and FS from the experimental isotope shift. The nonrelativistic NMS has been 
evaluated as the — m(l/Mi — I/M2) fraction of the nonrelativistic transition energy, obtained 
as the difference between the experimental value of the transition energy (16956.2 cm~^) and 
the related relativistic correction (45.4 cm~^), calculated in this work. As a result, the nonrel- 
ativistic NMS was obtained to be -547.5 MHz. The sum of this value and the relativistic NMS 
correction (-9.7 MHz), calculated in this work, yields -557.2 MHz for the total NMS value. 
Subtracting this NMS value and our theoretical FS value (6.9 MHz) from the experimental 
isotope shift, we obtain the value -208.2 MHz for the SMS and -105.8 (GHz amu) for the SMS 
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contribution to the constant K. This value differs from our theoretical result -98.5 (GHz amu) 
(Table 1) by 7% and from the result of [15] by 9%. The disagreement between the theoretical 
and experimental values is mainly caused by the correlation effects, which are extremely large 
for neutral Na. It is very difficult to compute the isotope shift parameters for neutral Na with 
high accuracy. For systematic investigations of the nonrelativistic isotope shift by the multi- 
configuration Hartree-Fock method we refer to [19], where an extremly large CI expansion was 
used. Although in our relativistic calculations the CI expansion is not so large, it is sufficient 
to show the role of the relativistic effects in the series of Na-like ions. Adding our relativistic 
correction to the nonrelativistic contribution from [19] yields K — —103.8 (GHz amu) that is 
in a fair agreement with the experimental value presented above. 

Table 2 shows the field-shift constant F, defined by SE-ps = —F6{r'^), for the 3pi/2 — 3s 
transition in Na-like ions. The comparison of the nonrelativistic (second column) and relativistic 
(third column) data demonstrates the role of the relativistic effects in the field isotope shifts. 
Our relativistic field-shift constants are in a good agreement with the related data from [15], 
presented in the fourth column of the table. 

In Table 3 we present the results of our calculations for the isotope shifts in 40^^,13+^36^^.13+^ 
4o^j,i4+^36^j,i4+^ and 86j^j,22+^84j^j,22+^ calculatc the field isotope shift, we used the rms nu- 
clear charge radii given in [20]. Since we consider transitions between levels which differ from 
each other only by the total angular momentum, the corresponding transition energies are com- 
pletely determined by relativistic and QED effects. In particular, it means that all the mass and 
field shift contributions given in the table are of pure relativistic origin. Comparing the individ- 
ual mass shift contributions given in the Table 3, we again observe a very significant cancellation 
of the contributions from the nonrelativistic part of the recoil operator (-f^^MS + -^sms) 
corresponding contributions from the relativistic part of the recoil operator {H^ug + -fastis)- 
can also be seen that the mass isotope shift strongly dominates over the field isotope shift for 
the ions under consideration. 

In Table 4, we compare the isotope shifts calculated in this paper with the total theoretical 
and experimental values for the transition energies under consideration [1]. The non-QED 
parts of the theoretical values for the transition energies were obtained by the large-scale CIDF 
method described in the previous section. The QED corrections were evaluated by using the 
one-electron Lamb shift data taken from [21] with an effective nuclear charge number Z^s- For 
a given one-electron state, Zes was chosen to reproduce the related DF electron charged density 
at the Compton wavelength distance from the nucleus. For the argon ions, these results are 
also in a fair agreement with the corresponding results of Refs. [22, 23] presented in the third 
column. In the sixth cohimn of Table 4 we present the total isotope shifts in 40^^,13+^36 ^j,i3+_ 
These data were obtained with the Dirac multiconfiguration wave functions and their relative 
errors are determined by the relative errors in the wave functions. As one can see from the 
table, the isotope shift in 4o^j,i3+^36^j.i3+ jg order of magnitude larger than the current 

experimental error. This provides very good perspectives for a first test of the relativistic theory 
for the recoil effect in middle- Z ions. 
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Table 1: The specific- mass-shift contribution to the isotope shift constant K (GHz amu) 
the 3pi/2 — 3s transition in Na-hke ions. 



Ion 


NSMS 


SMS 


SMS - NSMS 


RSMS 


Total 


SMS [15] 


Na 


-100.1 


-101.0 


-0.9 


2.5 


-98.5 


-97 


Mgi+ 


-412.7 


-416.1 


-3.4 


10.0 


-406.1 


-362 


A12+ 


-926.3 


-934.8 


-8.5 


24.0 


-910.8 


-837 


C16+ 


-5082.9 


-5162.4 


-79.5 


183.2 


-4979.3 


-4846 


Ar^+ 


-6648.2 


-6765.7 


-117.5 


261.6 


-6504.1 


-6363 


Pei5+ 


-26677 


-27713 


-1036 


1995.3 


-25718.6 


-25662 


Xe"3+ 


-200802 


-240996 


-40194 


69295 


-171701 





NSMS - nonrelativistic specific-mass shift 

SMS - specific-mass shift calculated with the Dirac wave function 
RSMS - relativistic-operator correction to the specific mass shift 
Total - SMS+RSMS 
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Table 2: Field shift constant F (MHz/fm^) for the 3pi/2 — 3s transition in Na-like ions. 
(r2)i/2 = 0.836^1/3 ^ 0.570 {A is the atomic mass number). 



Ion 


NFS 


FS 


FS [15] 


Na 


34.6 


36.45 


38.4 




115.9 


123.2 


125.8 


A12+ 


247.3 


265.3 


268.4 




1538.8 


1726.2 


1729.7 


Ar^+ 


2133.8 


2423.4 


2426.5 




14095.9 


18107.9 


18115 


Xe^3+ 


373760 


965738 





NFS - nonrelativistic field shift constant 
FS - relativistic field shift constant 
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Table 3: Individual contributions to the isotope shifts of the forbidden transition in ^°Ar/^^Ar 
and ^^Kr/^^Kr (cm~^). The rms nuclear charge radii used in the calculation are (r^)^/^ = 3.390, 
3.427, 4.188, and 4.184 fm for ^^Ar, ^°Ar, ^^Kr, and ^^Kr, respectively [20]. 



Ion 


Transition 


NMS 


SMS 


RNMS 


RSMS 


FS 


Total 


Ari3+ 
Kr22+ 


2s22pl 2p^/2 P3/2 

2si2pi ^Pi -3 P2 

3s23p2 3p, -3 P2 


0.1053 
0.0797 
0.0053 


-0.0742 
-0.0698 
0.0010 


-0.0822 
-0.0627 
-0.0025 


0.1151 
0.0887 
0.0000 


-0.0005 
-0.0001 
0.0001 


0.0635 
0.0358 
0.0039 



NMS - normal mass shift calculated with the Dirac wave function 
SMS - specific mass shift calculated with the Dirac wave function 
RNMS - relativistic-operator correction to the normal mass shift 
RSMS - relativistic-operator correction to the specific mass shift 
FS - field shift 
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Table 4: Energies of forbidden transitions for Ar and Kr ions and the isotope shifts in ^°Ar/^^Ar 
and ^^Kr/s^Kr. 



Ion 


Transition 


Theory [22, 23] 


Theory [1], this work 


Experiment [1] 


Isotope 






A (nm, air) 


A (nm, air) 


A (nm, air) 


shifts (nm) 




2s22pi 2p,/2 - 'P3/2 


440.99 


441.16(27) 


441.2559(1) 


0.00126 


^j,14+ 




593.88 


594.24(30) 


594.3880(3) 


0.00136 


Kr22+ 






383.35(95) 


384.1146(2) 


0.00005 
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